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Using the strong coupling expansion we calculate the analytically exact order parameter for the effective theory of 
Polyakov loops with applied chemical potential in the large- A"c limit. 
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1. Introduction 

In recent years the order parameters of Quantum Chro- 
modynamics at finite temperature and chemical potential 
(/i) have gained renewed attention. At least part of this 
comes from the need to know where and when phase tran- 
sitions occur in quark matter experiments, such as the 
observations of quark-gluon plasma at RHIC and LHC- 
ALICE. 

However, due to the confining nature of QCD, the usual 
perturbative approach breaks down at low energies, mean- 
ing that other methods need to be employed for quantative 
results in this regime. A widely succesful method is the 
lattice formulation of QCD [1 proposed by Wilson in 1974 
in which the challenges lie in the evaluation of group inte- 
grals over SU{i). 

Aside from the significant challenges in the analytical 
evaluation of these integrals, numerical calculations also 
run into trouble in the form of the sign problem when the 
chemical potential is applied. As such it is instructive to 
work in the — >■ oo limit, which, due to the \/Nc ex- 
pansion of QCD, provides a reasonably accurate approxi- 
mation to Nc = 3 but allows for a simplification of some 
problems. Here we will present analytically exact calcula- 
tions within the effective theory of Polyakov loops in the 
Nc — oo limit and determine an expression for the order 
parameter of QCD in the case of an applied real quark 
chemical potential. 

2. Lattice Formulation 

Using the lattice formulation for a pure gauge theory 
proposed by Wilson fp the Yang-Mills action can be writ- 
ten as 



A. 

2Nc 



^Tr[f/p + C/t] 



(1) 



where /? 



2Na 



and Up — UijUjkUkiUii are the plaque- 



tte variables taking their values in the fundamental rep- 
resentation of SU{Nc). Quarks are introduced by using 
Wilson's method of heavy fermions by adding (in the no- 
tation by Creutz [5]) 



where the sum is over all nearest-neighbour pairs 

with one term for each ordering of i and j, is a unit 
vector from i to j and Kf = g^J^^.^l_^_]J is the flavour- 
dependent hopping parameter with D being the number of 
spatial dimensions. At this point it is straightforward to 
integrate out the fermionic degrees of freedom and expand 
the fermion determinant in terms of Kf 



det(l - KfM{U)) exp - 



OO 



Tr[M^(J7)] 



(3) 



where the sum is over all loops without backtracking. A 
chemical potential is added by multiplying the positively 
(negatively) oriented time-like products by e""^ (e~°'^) in 
likeness of a fourth component of the gauge field [3]. At 
this point the effective theory of Polyakov loops can be 
derived using the strong coupling expansion together with 
the option of gauge fixing any maximal tree of link vari- 
ables Uij = 1 [1] or by using a Kadanoff-inspired bond- 
moving procedure [5]. To leading order in both (3 and n 
the theory takes the form 

S'e//(J,/ie°^-^,/ie°^-'') 

+ W\x)W{x + j)] 
+7V,^[/ie'^^-^iy(x) + her-'^-^'W^x)] 



(4) 



where W{x) = Tr J| U^i is the thermal Polyakov loop, 

n=0 

J > runs over the nearest neighbours only and J = 
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2{P/2N^)^- and h = 2^k^- to leading order for Nj 
mass degenerate quark flavours. Here is the number of 
lattice sites in the temporal direction. The values of J and 
h get renormalized by taking into account the spatial link 
integrations. For notational simplicity we define /i = aiV^/i 
from now on. The purpose of this paper is to derive the 
analytically exact expectation value, (W), of the Polyakov 
loop in the Nc — >■ oo limit, keeping the ratio ^ fixed. 



3. The Expectation Value 

The discussion below is largely analoguous to that of 
[6] but with a few crucial exceptions in order to treat the 
case of an applied chemical potential. Starting from the 
action Q and writing 



with 



S{J,0,0)^IjJ2{{W{x)}W\x + j) 



+ W{x) {W^ (x + j))- {W{x)) {W^ (x + j)) (5) 
+ [W(x) - {W{x))][W\x + j) - {W^x + j))] + h.c} 

the last term does not contribute in the Nc — > oo limit due 
to factorization. Translational invariance allows the j-sum 
to be performed for arbitrary number of spatial dimensions 

S{J,he^M'^) - NcYyhe^'W{x)^he-^W\x)] (6) 
+JDYy{W)W\x) + {W^)W{x) - {W){W^)] 

X 

Notice here that {WY ^ (W^^) since the action is non- 
Hermitian due to the chemical potential being real. In 
fact = {W) and (Pl/t)t = (p^/t) due to the Haar 

measure, meaning both expectations are real. The action 
can also be expressed as 

S{J,he^\he-^') = S{0,h+,h_) - JD^{W){W'^) (7) 



with h± given by the self-consistency relations: 

h-='^{W)+he-^ 

As usual, the expectation values are given by 
1 



(8) 
(9) 



{W) = 



z 



dW {We^p[S{0,h+,h^)]) (10) 



SU(N) 



1 

z 



exp[S{0,h+,h^)]) (11) 



SU(N) 



Z{J,h+,h^) ^ N{J) J dWeMS{0,h+,h^)] (12) 

SU{N) 



= N{J) J dWexp 



SU{N) 



where N{J) = eiip{-JDJ2A^){W^)) is a J-dependent 
prefactor. 

4. SU(7Vc) to U(iVe) 

In the large- A'c limit SU(A'c) integration is equivalent 
to U(-/Vc) integration. The relevant integral, as calculated 
by H, is 



det[Ar'/,-i(AO] 



U{N) 

rN-1 ■ 



A(A2) 



(13) 



where A? is the eigenvalues of AB and A(a::) is the Van- 
dermonde determinant 



N 



(14) 



In the relevant case A ^ h+I and B = h^I where / is an 
Nc X Nc identity matrix whi ch m eans that both nomina- 
tor and denominator of eq. (flSl have ^N{N — 1) zeros. 
Removing these using L'HopitaPs rule reduces the integral 
to 

/(/i+,/i_,0) = J dW exp[Nch+W + Nch^W''] (15) 

C/(Wc) 

= det[Ij^,{2Nc^h+h^)] (16) 

up to an overall, and in our case irrelevant, normalization 
constant. Amazingly enough this equation means that the 
integral, /, only depends on the product of h+ and h- , sim- 
plifying the problem of calculating the expectation values 
markedly. Reexpressing the partition function: 



Z{J, h+, h^) = Z{J, ^h+h^, V^+h-) (17) 
means that the expectation values can be expressed as 
1 d 



(W) 



NcN^Z dh 

1 [hZ d 
2NcN^.Z\j h^dg 
1 d 



Z{J,9,g) (18) 



NcN^Z dh 
2NcN^Zy ~hZdg 



Z{J, y/h+h_, y/h+h-) 

Z{J,9,g) (19) 
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for g = yJhZh^^ 7^ and = The case 

of = and = is trivially solved and yields 
(W) — (W^) = 0. Notice that from these expressions 
h+{W) = h-{W^), which, by using the self-consistency 
relations, also leads to h^he~^ 
integration can be performed: 



/i_/ie''. Using 015] the 



h+{W) _ h^{W^) 



h^h^ for yJTuJi- 



Nr. 



Nr. 



for y/h^h_ 



< 

> A 



(20) 



Solving the self-consistency equations in the case of yJhj^KZ 
< i is straightforward, but in the case of yjhj^h^ > ^ 
there is a small detail which requires mentioning. Rewrit- 



ing the integration result ( 20 ) 



(W) 





(21) 



we can use that for h ^ then ^ = e^^ to simplify the 
expression. To make this step cover the case of ft. = 
requires setting fi = simultaneously, such that ^ = 1 
which is the solution when h = from the beginning. 
Now, having taken the preliminary steps, the rest of the 
derivation is simple and leads to the solutions 



Nr 

and 

Nr 




exp[-A] 



(22) 




exp[A] 



(23) 



with a transition point between the two solutions at JD = 
1 — 2h. As can be seen from the solutions, simply setting 
h = does not remove the /i-dependence; it is necessary 
to also set n = 0, revealing that h — cannot be reached 
in the /i — > limit of the solution. Figure [l] showing the 
convergence of the numerically calculated (W) values for 
U {Nc) when subjected to the self-consistency relations in 
the case of A'c = 1, 3, 7, 15 together with the analytical ex- 
pression in the Nc ^ 00 limit. The associated free energy 
is found to be 



F(J) 



JD{W){W^)-h+h^ 
JD{W) (M^t) _ + \ ln(2yMi^ 



(24) 



with a third order phase-transition point at JD = 1 — 2/i in 
which {W) — ^ exp(/i) and (W^) = ^exp(— /i). The tran- 
sition line is illustrated in fig. [2] Inserting the solutions 
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Figure 1: Plot of (W) /Nc for h = 0.2 and jl = 0.4 for Nc = 1,3,7, 15 
and 00 in the case of U{Nc). 



Free Energy 




Figure 2; Plot of the free energy F(JD,h) with the transition line at 
JD = l-2h added. 



the JD < 1 ~ 2h part of the free energy can be simplified 
and written as 



F{J) = 



1 - JD 



for JD <l-2h 



(25) 



The free energy, like the partition function in |17[ is inde- 
pendent of the chemical potential, since the factors of e'* 
and e^^ cancel - a consequence of this is that the quark 
number density is zero 



d_ 



(26) 



This result is an artifact of U{Nc) and not SU{Nc) - a 
claim which is supported in a recent mean field study j9j 
further investigating the U{Nc) and SU{Nc) difference. 
Note that, surprisingly, this means that the N^ — > 00 limit 
of the U{Nc) and SU{Nc) theories differ even though the 
difference between the two groups vanishes as l/N^- 

For recent numerical work for SU(3) in the case of a 
real chemical potential, we refer the reader to [101 HB IE] 
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where complex Langevin techniques and the flux repre- 
sentation is used. Work on a related model incorperating 
higher order contributions to the effective theory is found 
in [13] and [13] respectively. It is remarkable that the 
phase diagram of the Nc = 3 theory [14] looks qualitively 
similar to that of fig. [2] 

5. Conclusion 

We have found an analytically exact expression for the 
order parameter of QCD, {W), within the effective theory 
of Polyakov loops with an applied chemical potential in 
the limit Nc oo. Furthermore the free energy has been 
calculated and displayed a third order phase transition at 
JD = 1 — 2h. The free energy, and hence also the point 
of the phase transition, was found to be /t independent, in 
agreement with large- iVc counting. The expectation values 
of the Polyakov loops display the dependence on quark 
chemical potential that could be expected from intuitive 
arguments. 
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